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Dark energy (i.e., a cosmological constant) leads, in the Newtonian approximation, to a repulsive 
force which grows linearly with distance and which can have astrophysical consequences. For exam¬ 
ple, the dark energy force overcomes the gravitational attraction from an isolated object (e.g., dwarf 
galaxy) of mass 10 7 M@ at a distance of 23 kpc. Observable velocities of bound satellites (rotation 
curves) could be significantly affected, and therefore used to measure or constrain the dark energy 
density. Here, isolated means that the gravitational effect of large nearby galaxies (specifically, of 
their dark matter halos) is negligible; examples of isolated dwarf galaxies include Antlia or DDO 
190. 


I. INTRODUCTION 

The discovery of dark energy, which accounts for the 
majority of the energy in the universe, is one of the most 
significant of the last 20 years. While the repulsive prop¬ 
erties of dark energy are well known in the cosmolog¬ 
ical context, they have not been as thoroughly under¬ 
stood on shorter, astrophysical, length scales. Previous 
work has constrained the cosmological constant on solar- 
system scales [I], but its effects are obviously too small 
to be directly observed. 

In what follows, we discuss the repulsive dark energy 
force and its astrophysical effects on galactic scales. Be¬ 
cause this force grows linearly with distance, its effect is 
most likely to be significant for weakly-bound satellites 
with large orbits. To detect the effects of dark energy, we 
must first understand orbits resulting from ordinary grav¬ 
itational dynamics, with potentials mainly determined by 
the distribution of dark matter. We find that observa¬ 
tions of distant satellites of isolated dwarf galaxies could 
be used to detect the effects of dark energy. Here, iso¬ 
lated means sufficiently far from other sources of gravita¬ 
tional potential. When dwarf galaxy systems are not suf¬ 
ficiently isolated, the orbits of their satellites are subject 
to tidal forces from nearby large galaxies. These tidal 
forces can distort orbital shapes, and enforce an upper 
limit on orbital radii. 


II. NEWTONIAN GRAVITY AND 
COSMOLOGICAL CONSTANT 

The Einstein equation with cosmological constant A is 

Rpv - 2 9pvR = 8-itGT^ + g /J ,„ A. (1) 

Contracting both sides with g M1/ , one gets R = —8 7r G T— 
4 A where T = is the trace of the matter (including 


dark matter) energy-momentum tensor. This can be sub¬ 
stituted in the original equation to obtain 

Rp V = 8 7T G ^ i g^ v T ^ - g^ A . (2) 

In the Newtonian limit, one can decompose the metric 
tensor as g^ + h^ with \h^„\ -C 1. Specifically, 

we are interested in the 00-component of the Einstein 
equation. We parameterize the 00 th -component of the 
metric tensor as 

ffoo = 1 + 2 $ , (3) 

where $ is the Newtonian gravitational potential. To 
leading order, one can show that Q 

Roo ~ 2 5°° = V 2 $ . (4) 

In the inertial frame of a perfect fluid, its 4-velocity is 
given by = (1, 0) and we have 

Tpv = (p + p)u tl u v - pg^ = diag( p,p), (5) 

where p is the energy density and p is the pressure. 
For a Newtonian (non-relativistic) fluid, the pressure is 
negligible compared to the energy density, and hence 
T k Too = P- As a result, in the Newtonian limit, the 
00-conrponent of the Einstein equation reduces to 

V 2 $ = 4ttGp-A, (6) 

which is just the modified Poisson equation for Newto¬ 
nian gravity, including cosmological constant. This equa¬ 
tion can also be derived from the Poisson equation of 
Newtonian gravity, V 2 $ = 47rG(p + 3 p), with source 
terms from matter and dark energy; p « 0 for non- 
relativistic matter, and p = —p for a cosmological con¬ 
stant. 

Assuming spherical symmetry, we have V 2 $ = 

r? ~§r ( r2 ) an d the P°i sson equation is easily solved 
to obtain 
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( 7 ) 
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where M is the total mass enclosed by the volume | 7 tr 3 . 
The corresponding gravitational field strength is given by 


9 = -V$= 



r . 


( 8 ) 


Therefore, the cosmological constant leads to a repulsive 
force whose strength grows linearly with r. 

One can also derive <7 by starting with the de Sitter- 
Schwarzschild metric 


ds 2 = 


1 - 


2 GM A 


- —r dt 


( 9 ) 
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2 GM A 2 

-T r 

r 3 
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dr 2 - r 2 dil 2 


Galaxy Mass r c 


1O 6 A7 0 

10.7 kpc 

io 7 al 0 

23.1 kpc 

io 8 ai 0 

49.8 kpc 

io 9 al 0 

107 kpc 

lO lo Af 0 

231 kpc 

10 11 Mq 

498 kpc 

io 12 ai 0 

1.07 Mpc 

io 13 m 0 

2.31 Mpc 

io 14 ai 0 

4.98 Mpc 


TABLE I: Galaxy masses (units of solar mass Mq) and the 
corresponding r c . 


which describes the spacetime outside a spherically sym¬ 
metric mass distribution M in the presence of a cosmo¬ 
logical constant A. One then obtains Eq. 0 and hence 
Eq. (JSI) by identifying Eq. (J3J) with the 00-component of 
the de Sitter-Schwarzschild metric. 


III. GALAXIES 

The results obtained in the previous section are rele¬ 
vant to galaxies. For instance, in the presence of the cos¬ 
mological constant A, Eq. ([ 8 ]) describes the Newtonian 
gravitational field strength outside a galaxy with a spher¬ 
ically symmetric mass distribution M. From Eq. ([ 8 |), it 
is clear that when r is sufficiently large, the repulsive 
dark force will dominate over the gravitational attrac¬ 
tion. The critical value of r beyond which this happens 
is given by 


/ 3GM \ 1/3 _ / 3 M \ 1/3 

V A ) [sttpa) 


( 10 ) 


where p\ = ~ (2.3 x 10 _3 eV ) 4 is the observed 

energy density of the cosmological constant. Table I dis¬ 
plays galactic masses in units of solar mass Mq and the 
corresponding r c . 

Typical galaxies, including our Milky Way, have total 
mass (including dark matter) > 1 O 11 _ 12 M 0 and sizes 
~ 50kpc. According to Table I, r c > 500 kpc for these 
galaxies, so the dark force is not likely to affect internal 
dynamics, but may impact galaxy-galaxy interactions ij], 
and limit the size of galaxy clusters (~ 10 14 Mq, size ~ 
Mpc). 

Some dwarf galaxies have total mass (including dark 
matter) ~ 10 7 Mq. These include Ursa Major II, Coma 
Berenices, Leo T, Leo IV, Canes Venatici I, Canes Ve- 
natici II, and Hercules (analyzed by 1 ), and also Leo II 
(1 and Leo V @]. The irregular galaxies Leo A [§}, Antlia 
| and DDO 190 0 also have masses around 10 7 Mq. 

For galaxies with mass ~ 10 7 Mq, we have r c ~ 23 kpc. 
Thus, their galactic rotation curves could be affected by 


the dark force: rotational velocities of stars or gas clouds 
bound to these galaxies should be smaller than that pre¬ 
dicted by ordinary Newtonian gravity. This in turn could 
provide a novel way to measure the cosmological con¬ 
stant in the future. Rotation curves for many galaxies 
have been measured to radii of ~ 30 kpc or more, and for 
some dwarf galaxies to ~ 10 kpc Q Low surface bright¬ 
ness (LSB) galaxies may also be worthy of investigation 
[jjj]- Some LSBs with total mass ~ 10 lo Mg have disks 
as large as 100 kpc. 

The Navarro-Frenk-White (NFW) profile [l3[ is a com¬ 
monly used parametrization of dark matter halo energy 
density: 


r/Rs{l + r/R s f' 

where p 0 is a characteristic halo density and R s is the 
scale radius. These two quantities vary from galaxy to 
galaxy. While the detailed shape of the actual dark mat¬ 
ter density may differ from the NFW profile, the asymp¬ 
totic 1/r 3 behavior is widely accepted. Our results below 
will not be sensitive to the density profile at small r. 

Consider a dwarf galaxy (DG) and a larger galaxy (LG) 
(e.g., the Milky Way) whose centers of mass are separated 
by a distance R , and a satellite of the DG whose orbital 
radius is roughly r. If the distance R is sufficiently large, 
we can neglect the gravitational potential of the LG and 
treat the DG-satellite system as approximately isolated. 
In that case, the values in Table I provide a rough guide 
for distances r at which the dark energy force becomes 
significant. In the following section we will investigate 
to what extent measurement of satellite velocities can 
constrain the dark energy density around the DG. 

But first let us examine in more detail under what cir¬ 
cumstances we can neglect the gravitational effects from 
the (dark matter halos) of neighboring galaxies on the 
DG. We will assume an NFW profile for both the DG 
halo and the larger galactic halo. The distance R from 
the center of the DG to the center of the LG is generally 
not equal to the distance from the satellite to the cen¬ 
ter of the LG, which can vary from (R — r) to (R + r). 
Therefore, the gravitational pull exerted on the satellite 
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Dwarf Galaxies distance to MW distance to M31 Dwarf Galaxies negligible tidal effect 


Leo T 

422 kpc 

991 kpc 

Leo T 

r 

< 

5-6±l;| kpc 

Leo IV 

155 kpc 

899 kpc 

Leo IV 

r 

< 

l-9tj ;l kpc 

Canes Venatici I 

218 kpc 

864 kpc 

Canes Venatici I 

r 

< 

2 - 6 ii j kpc 

Canes Venatici II 

161 kpc 

837 kpc 

Canes Venatici II 

r 

< 

2.2^1 kpc 

Hercules 

126 kpc 

826 kpc 

Hercules 

r 

< 

L4±o.l k P c 

Leo II 

236 kpc 

901 kpc 

Leo II 

r 

< 

3-lths kpc 

Leo V 

179 kpc 

915 kpc 

Leo V 

r 

< 

2 - 8 ±i;f kpc 

Leo A 

803 kpc 

1200 kpc 

Leo A 

r 

< 

10 . 6 tg 2 g 6 kpc 

Antlia 

1350 kpc 

2039 kpc 

Antlia 

r 

< 

IS-Stme k P c 

DDO 190 

2793 kpc 

2917 kpc 

DDO 190 

r 

< 

39.4 tlH kpc 


TABLE II: Some dwarf galaxies with mass ~ 10 7 Mq and 
their distances from the Milky Way (MW) and the An¬ 
dromeda galaxy (M31). See Table 2 in [l5) |. 


Dwarf Galaxies 

^half 

po / GeV cm 3 

Leo T 

178 ± 39 pc 

{0.028,0.22, 1.74} 

Leo IV 

116 ± 30 pc 

{0.037,0.25,1.66} 

Canes Venatici I 

564 ± 36 pc 

{0.017,0.21,2.69} 

Canes Venatici II 

74 ± 12 pc 

{0.053,0.30,1.65} 

Hercules 

330 ± 63 pc 

{0.020,0.20,2.05} 

Leo II 

123 ± 27 pc 

{0.035,0.24,1.67} 

Leo V 

42 ± 5 pc 

{0.086,0.38,1.72} 

Leo A 

354 ± 19 pc 

{0.019,0.20,2.11} 

Antlia 

471 ± 52 pc 

{0.018,0.21,2.42} 

DDO 190 

520 ± 49 pc 

{0.017,0.21,2.56} 


TABLE III: Some dwarf galaxies with mass ~ 10' Mg and 
their half-light radii rhaif and po- For each galaxy, the three 
different values of pa correspond to the exponents {— 1.6 — 
0.4, -1.6, -1.6 + 0.4} in Eq. (flG]). 


by the LG is different from the pull on the DG, leading 
to a tidal effect. (See 0) for previous work regarding 
the tidal effects on orbiting satellites around their host 
galaxies.) This tidal effect is repulsive: it pulls apart the 
DG-satellite system. Perhaps surprisingly, for many DGs 
(i.e., near the Milky Way), the tidal effect is large enough 
to distort and even destabilize the satellite orbits. 

Let the total mass of the DG enclosed within r be 
-L/dg(V) and that of the LG enclosed within R ± r be 
± r). Then we have 


M DG (r) = / Airr' 2 p DG (r')dr', (12) 

Jo 

/*i?±r 

M LG (R±r) = / Anr' 2 p LG (r')dr '. (13) 

Jo 

The circularity and stability of the satellite orbits can be 
guaranteed by requiring that the magnitude of the tidal 
force due to the LG, F^q, is much smaller than the 


TABLE IV: Required orbital radii for the satellites of some 
dwarf galaxies with mass ~ 1 O 7 M 0 to ensure negligible tidal 
effect. The upper limit on r is determined by requiring 
ptida.1 /f dq < o p The cen tral value corresponds to the expo¬ 
nent — 1.6 in Eq. © while the ± values correspond to the 
exponents —1.6 ± 0.4. 


gravitational pull due to DG, F D q. This requires 


TTitidal 

^LG 


G Mlg(R) r 
R 2 R 


Fdg 


GM DG (r ) 

r 2 


, (14) 


which implies 


/ Mlg(-R) 
V M DG (r) 



« 1 . 


(15) 


According to 0 ]. many dwarf galaxies are at least 100 
kpc away from the Milky Way and much farther from the 
Andromeda galaxy (M31). Some of these dwarf galaxies 
with mass ~ 1O 7 M 0 include Leo T, Leo IV, Canes Ve- 
natici I, Canes Venatici II, Hercules, Leo II, Leo V, Leo 
A, Antlia and DDO 190. Their distances from the Milky 
Way and M31 are shown in Table II. 

For the Milky Way (MW), we have po ~ 0.2 GeV cm -3 
and R s ~ 25 kpc (see Fig. 1 in 0 )- For dwarf galaxies 
with mass ~ 10'Mg, the analysis in 0 suggests that 
R s ~ 0.795 kpc (see their Table 3 which gives the best fit 
parameters for some dwarf galaxies assuming the NFW 
profile). In terms of the half-light radius rh aif ( the radius 
at which half of the total light is emitted), [17j obtains a 
relation for the mean density {p) interior to rhaif: 

/ \ -1.6T0.4 

(p) ~ 2600 ( ^ ) GeV cm" 3 . (16) 

(Note that M 0 pc" 3 ~ 40 GeV cm" 3 ). In the second 
column of Table III, the half-light radii of some dwarf 
galaxies with mass ~ 1O 7 M 0 are listed. We adopt these 
values of rhaif from 00 0 | . For each of the galaxies 
listed in Table III, by using r = rhaif and the correspond¬ 
ing value of (p) obtained from Eq. (fl6l) (together with 
Rs ~ 0.795 kpc) in the NFW density profile, we estimate 
the value of po- The three different values of po listed in 
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N 

orbital radii 

95% Cl of c/10 -84 GeV 2 

5 

1-10 kpc 

{1.31, 1.78} 

5 

5-15 kpc 

{ 1.49, 1.66} 

5 

10-20 kpc 

{1.54, 1.61} 

5 

15-25 kpc 

{ 1.55, 1.60} 

5 

20-30 kpc 

{ 1.56, 1.59} 

10 

1-10 kpc 

{ 1.43, 1.70} 

10 

5-15 kpc 

{ 1.53, 1.64} 

10 

10-20 kpc 

{1.56, 1.61} 

10 

15-25 kpc 

{ 1.56, 1.60} 

10 

20-30 kpc 

{ 1.57, 1.59} 

TABLE V: 
error in a 2 . 

95% confidence interval (Cl) of c, assuming 
N is number of satellites. 

N 

orbital radii 

95% Cl of c/10 -84 GeV 2 

5 

1-10 kpc 

{0.40, 2.31} 

5 

5-15 kpc 

{1.11, 1.81} 

10 

1-10 kpc 

{0.68, 2.07} 

10 

5-15 kpc 

{ 1.25, 1.78} 


TABLE VI: 95% confidence interval (Cl) of c, assuming 5% 
error in v 2 . N is number of satellites. 


the third column of Table III correspond to the exponents 
{-1.6 - 0.4, -1.6, -1.6 + 0.4} in Eq. (HU). 

Using the normalization factors (po and R s ) for both 
the Milky Way and the dwarf galaxies with mass ~ 
1O 7 M 0 , we can determine the satellite radii for which 
Eq. (fl5l) is satisfied. In Table IV, we display the re¬ 
quired orbital radii for the satellites of some dwarf galax¬ 
ies with mass ~ 10'M 0 , assuming negligible tidal ef¬ 
fect. The upper limit on r is determined by requiring 
-Flq^/Fdg <0.1. The central value corresponds to the 
exponent —1.6 in Eq. (USD while the ± values correspond 
to the exponents —1.6 ± 0.4. 


IV. CONSTRAINTS ON DARK ENERGY 
DENSITY FROM ROTATION CURVES 


For galaxies with mass ~ 1O 7 M 0 , we have r c ~ 23 kpc, 
so the effect of the dark energy force becomes observa- 
tionally significant for r ~ 20 kpc. For such radii, Table 
IV indicates that the tidal effect due to the Milky Way 
might have already distorted and even destabilized the 
satellite orbits in the dwarf galaxies Leo T, Leo IV, Canes 
Venatici I, Canes Venatici II, Hercules, Leo II, Leo V and 
Leo A. On the other hand, satellite orbits of Antlia and 
DDO 190 could be strongly affected by the dark energy 
force at radii for which the tidal effect due to the Milky 
Way is negligible. 

For isolated dwarf galaxies such as Antlia and DDO 
190, the rotational velocity-squared v 2 of their satellites 


at r is given by 


2 GMdgW 1 A 2 

v (r) =-— A r 


(17) 


r 3 

where MdgM can be obtained by a simple integration: 

r + R s 


M DG (r) = 4irp 0 R' i s 


In 


Rs 


r + R s 


■ ( 18 ) 


Thus, for a set of measurements on v 2 (r) at some level 
of sensitivity, we fit v 2 (r) with 


v 2 (r) = - 
r 


In 


r + b 


(19) 


where a and b are some constants. A (statistically sig¬ 
nificant) positive fit value for c suggests the existence of 
a cosmological constant. The cosmological value for c is 
c = 1.58 x 10 -84 GeV 2 . (Note that in a realistic situation 
we cannot rely on the NFW parametrization being exact, 
and the fitting function should probably be slightly more 
general than the one used above.) 

In Table V, we simulate the results of measurements 
on v 2 (r ) with corresponding error of 1%. We take 
Po ~ 0.2 GeV cm -3 and R s ~ 0.795 kpc for the dwarf 
galaxies. We vary the number of satellites N and their 
(randomly generated) orbital radii. For example, at 95% 
confidence level, one could bound c to be positive using 
5 satellites at r ~ 1 — 10 kpc. In order to bound c close 
to its cosmological value, one would need, e.g., at least 5 
satellites at r ~ 10 — 20 kpc or 10 satellites at r ~ 5 — 15 
kpc. 

In the event that dark energy is dynamical [l9j ]. as 
opposed to a rigid cosmological constant, it might form 
inhomogeneous clumps on galactic length scales. This 
behavior could, in principle, be detectable through the ef¬ 
fects discussed here: A from astrophysical measurements 
would be larger than the known cosmological value. In 
Table VI, we simulate the results from measurements 
on v 2 (r), assuming that the corresponding error is 5%. 
Again, we take po ~ 0.2 GeV cm” 3 and R s ~ 0.795 kpc 
for the dwarf galaxies. The table indicates that even 
at the sensitivity of 5%, one could rule out (at 95% 
confidence level) any A that is significantly larger than 
1.58 x 10" 84 GeV 2 by using, e.g., 5 satellites at r ~ 1 -10 
kpc. The very existence of satellites of dwarf galaxies 
(even those close to the Milky Way, and hence subject 
to significant tidal forces that limit r) provides an upper 
limit on the local dark energy density, probably no more 
than an order of magnitude larger than the cosmological 
value. 


V. MISSING SATELLITE PROBLEM 

Observations indicate fewer satellite galaxies than pre¬ 
dicted by numerical simulations involving cold dark mat¬ 
ter. This is known as the missing satellite problem [20j . 
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For instance, simulations predict a few hundred satellite 
galaxies within a few Mpc radius of the Local Group, but 
we have observed at least five times fewer. 

For the Local Group with mass ~ 10 13-14 M 0 , we have 
r c ~ a few Mpc, which naively suggests that the dark 
force might play a limiting role in the binding of satellite 
galaxies. However, if simulations of galaxy formation are 
performed with sufficient resolution and sufficiently late 
end-times in an expanding ACDM universe, the effect of 
the dark energy has already been accounted for (barring 
more exotic scenarios, like clumping of dark energy). 

In other words, we are not suggesting the dark energy 


force discussed here as a possible solution to the missing 
satellite problem. We are instead suggesting a possible 
alternative method for measuring the local dark energy 
density through rotation curves of dwarf galaxies. 
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